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Abstract 



The recent Super-Kamiokande data now admit only one type of mass hierarchy in a frame- 
work with three active and one sterile neutrinos. We show that neutrino masses and mixings 
generated by R-parity- violating couplings, with values within their experimental upper lim- 
its, are capable of reproducing this hierarchy, explaining all neutrino data particularly after 
' including the LSND results. 

v3 ; PACS number(s): 12.60.Jv, 14.60.Pq, ILSO.Fs 

o 

I The vacuum oscillation interpretation of solar neutrino data requires Am^^jj^j. ~ 10^^'^ eV^, 

O ' while the matter enhanced MSW solution prefers the range Arri^^i^j. ~ (10~^ — 10~^) eV^ 

. The atmospheric neutrino anomaly can be explained by Am'^^^ ~ (2 — 5) x 10~^ eV^ (with 

i-C I sin^ 2^atm > 0.88) [@, ^. In the standard three neutrino framework, which offers two independent 

mass differences, it is possible to conceive of mass hierarchies which can explain both solar and 
Oh. atmospheric results. If in addition the LSND data, carrying a positive indication of <-> 

(— I I oscillation with Aml^^Y^ ~ (0.3 — 1) eV^ are also sought to be simultaneously explained, 

then one has to expand the standard three neutrino scenario as it falls short of one independent 
mass difference. In fact, within this three neutrino framework one can fit any two of the three Am^ 
mentioned above. Under the circumstances, one often leaves the LSND data out of consideration 
' pending further confirmation from the MiniBooNE at FNAL or MINOS long baseline |p 

experiments, since KARMEN ^ does neither confirm nor exclude the LSND results. On 
the other hand, if one includes the LSND results to be explained together with the solar and 
atmospheric data, the minimal extension of the standard scenario that needs to be done is to 



add one sterile neutrino (ug) [^, 10 1 to the list of three active states. The resulting four neutrino 
picture can explain all the data [pd| ]. 

What are the possible choices of mass hierarchies among these four neutrinos? The data 
suggest that the choices are very limited. First, we devide them into two types, I and II. In type 
I, there are three almost degenerate states explaining the solar and atmospheric results, with a 
fourth state separated by a large gap. This hypothesis cannot explain the LSND results, since 
the LSND data require that the separated fourth state will have to be either z^e or v^, but each 
of them will have to be closely spaced with a third state to explain the solar and atmospheric 
results. If one disregards the LSND results till further confirmation, then of course type I scenario 
is allowed with the isolated state either I's or Vr- In type II, there are two pairs of approximately 
degenerate states separated by a large gap. The mixing between the two pairs is very small. 
This scenario can explain all the data, or in other words, the large gap could be the LSND 
gap. Two cases may arise in this framework. In type Ila scenario, (z^t— i^g) form one pair which 
explains the solar neutrino data, while {v^-Vs) form the other pair maximally mixed to explain the 
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atmospheric anomaly. In type lib, (i^r-'^fj.) form the pair that explains the atmospheric anomaly, 
while (us-t^e) pair explains the solar data Q. The most recent Super-Kamiokande (SK) data rule 
out the ^ Us interpretation of atmospheric neutrino anomaly at 99% CL j^], thus strongly 
disfavouring the type Ila scenario. We are then left with type lib as the only surviving option, 
shown in Fig. 1^. It should be noted that the oscillation data cannot discriminate between cases 
that occur by interchanging either the members within a given pair or the relative spacing of the 
two pairs. These cases are not separately shown. 

Under the circumstances, it is a timely exercise to identify those scenarios which reproduce 
the type lib spectrum. Do the R-parity- violating (^) supersymmetric models [14] fall in that cat- 
egory? In this note we seek to find an answer to this question. Defined in terms of lepton number 
(-L), baryon number (B) and spin (S) of the particle as i? = (— l)3-B+L+25^ -j-j^jg (^jgcj-g^g parity 
is +1 for all SM particles and -1 for their superpartners. Since neither L- nor ^-conservation is 
ensured by gauge invariance or any such fundamental principles, R-parity is an ad hoc symmetry 
put in by hand. Although there is no experimental confirmation yet in favour of non-vanishing ^ 
interactions, the neutrino oscillation data are somehow suggestive that it would be premature to 
abandon those couplings, as the origin of neutrino masses and mixings could be traced to some 
of these non-vanishing L-violating couplings. In order not to allow rapid proton decay we do not 
switch on L and B violations simultaneously. The following L-violating terms in the superpo- 
tential are then allowed: (a) XijkLiLjE^, (b) X[jj^LiQjD'^, and (c) HiLiHu- In these expressions, 
Li and Qi are SU(2)-doublet lepton and quark superfields, Ef and Df are SU(2)-singlet charged 
lepton and down quark superfields, and is the Higgs superfield that gives masses to up-type 
quarks. Stringent experimental constraints on these a priori independent parameters exist in the 
literature ||l5|| . 

Attempts to fit the observed neutrino data by masses and mixing angles generated by ^ 
couplings have been done in the past, both in the context of bilinear |16| as well as trilinear 
17 1 L-violating parameters. Most of the analyses have been carried out in the three-neutrino 



picture. A few remarks on those analyses are in order. Most of them discarded the LSND data as 
something not so reliable, and confined the discussions to the possibility of fitting only solar and 
atmospheric data. Adhikari and Omanovic [17|, on the other hand, tried to fit solar, atmospheric 
[preliminary SK) and LSND all together in a three-neutrino picture and that too considering only 
trilinear ^ couplings. But they did not consider the SK zenith angle dependence and also assumed 
an energy independent solar neutrino solution by ignoring the Chlorine data of the Homestake 
mine experiment. In fact, if one takes all present data into consideration, it is not enough just 
to add bilinear ^ terms to their analysis - the data compel one to go for a four-neutrino picture. 
Ref. |l^] deals with the simultaneous presence of bilinear and trilinear parameters. Although the 
emphasis in ref. |18] is mostly on finding an explanation for the solar and atmospheric neutrino 
data in a three-neutrino framework, a qualitative discussion on how to simultaneously explain 
the LSND data by admitting a fourth sterile state has also been presented. 

In the present work, we perform a numerical study to examine whether R-parity violation, 
with bilinear and trilinear parameters together, can reproduce the type lib spectrum. Even 
though the data allow an interchange of the location of the two pairs, we work in a situation 
where {u^-Vt) form the heavier paii|^. We assume that the masses of all active neutrinos are 
generated by R-parity violation. Since the mixing between the two pairs is small, we can, for all 
practical purposes, focus on the heavier pair and work in the y^i-fr subspace, always assuming 

^U:^ ^ Us solar neutrino oscillation will be tested by the SNO experiment soon |l2[ . 

^This framework is consistent with the Big-Bang Nucleosynthesis constraints on sterile - active mixings [ p^ . 
•^This hierarchy is stable under radiative corrections [Q. 
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that the parameters A'^^ and fii, responsible for the Ve mass generation, are much smaller. We 
also assume that the sterile state receives mass of the order of the fg mass from some other 
source. More specifically, we concentrate on those parameters which generate the masses and 
mixing angles in the above (2 x 2) subsector, keeping in mind that each of the two absolute 
masses should be of the order of the LSND gap. In other words, we parametrize the v^-Vt mass 
matrix in terms of the 1^ couplings, vary them within physical ranges, and then observe whether 
there exist solutions that simultaneously satisfy the following experimental constraints: 

^"iisND ~ ml ~ mi ~ (0.3 - 1) eV^, 

Am2,^ = ml - mi ~ (2- 5) X 10-3 eV2, (1) 
sin^ 26iatm = sin^ 26*34 > 0.88. 

Now we turn our attention to the analytic expressions of the neutrino masses induced by 
bilinear and trilinear couplings. We write the mass matrix as My = 7W*'^^ + A4|,°°p. The bilinear 
couplings contribute to the tree mass |16|. In a basis where there are no sneutrino vacuum 
expectation values (VEVs) (2^, this can be expressed asQ 

-^l'" = 92- -l^il^i'v'd = afiifii' (vd = {m)) , (2) 

where Mi^2 are the gaugino masses, and det M is the determinant of the (4 x 4) neutralino 
mass matrix in the R-parity-conserving case. Considering only the A' couplings, the one-loop 
squark-mediated contribution to the neutrino mass, in the same basis as before, can be written 

m 



as 



/K/m|) ^ f{ml/ml) 



(3) 



where f{x) = (xlnx — 2; + l)/(x — 1)^. Here, rrid,. is the down quark mass of the ith generation, mj, 

is an average of du and dju squark masses, and A'c = 3 is the colour factor. While writing Eq. 
we assumed that the left-right squark mixing terms are family-diagonal and are proportional to 
the corresponding quark masses, i.e., rn1^^{i) = md^m^ . The expression of the A-induced slepton- 
mediated contributions to the neutrino mass is similar to Eq. (^), and we do not display it here^. 
In the basis (i/^, z^t), the mass matrix can be parametrized as 

_ / KXl + af4 KX2X3 + a/i2M3 \ (.-^ 
[ KX2X3 + 0^.2^3 KXl + ai^l )' 

where a is given by Eq. (|2|), A2,3 are two generic trilinear couplings, and K captures the loop 
factors that enter into Eq. (^). 

We observe that the determinant of the mass matrix with only bilinear or with only trilinear 
couplings is identically zero^ But this yields a big hierarchy between and m^, contrary to 

''ill principle, one can rotate away the bilinear fj,i terms from the superpotential, but stiU having the AL = 2 
effects via the presence of sneutrino VEVs after the minimization of the scalar potential. One can as well work in a 
framework where there are no sneutrino VEVs but the /li parameters are present in the theory [^o|. We work in the 
latter basis. We emphasize that going from one basis to another, the parametrization would change, but the general 
conclusion we draw at the end remains unaffected. For a discussion on the basis-independent parametrizations of 
R-parity violation, see refs. pH . 

^We neglect the one-loop diagrams induced by the product of bilinear and trilinear couplings p^ . 

^If we have more than one combination of trilinear couplings, then although the determinant of the mass matrix 
with trilinear couplings alone will not be zero, still much more fine-tuning may be necessary in order to arrive at 
any solution with only trilinear couplings. 
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the requirement that these two states should be approximately degenerate. Only after taking 
the bilinear and trilinear parameters together we obtain solutions that pass the test in Eq. @). 
Actually we need at least four independent input parameters in the mass matrix in order to fit 
the data. In this case, they are the two trilinear couplings (A2, A3) and the two bilinear mass 
parameters (/X2, /^s)- Notice that these four parameters as well as a can take either sign, while 
K is always positive. We point out that the tree and loop contributions may have the same 
or opposite signs, and this relative sign plays a crucial role in deciding which combinations of 
parameters are allowed by the data. After diagonalizing the mass matrix in Eq. (Q), we demand 
that the eigenvalues 7713^4 and the mixing angle 6*34 = ^atm satisfy Eq. (|l|). In Figs. 2a and 2b we 
have displayed only a part of the solutions by plotting the acceptable mass spectrum as a function 
of some allowed input parameters, just to demonstrate that the mechanism of R-parity violation 
works as a viable explanation]]. The above parametrization is rather general as one can apply it 
for any A' (or A, for that matter) couplings irrespective of the second and third generation indices. 
Also at this stage one need not specify the squark (or slepton) and gaugino masses as they are 
absorbed in K and a respectively. We observe that the 'filter' of Eq. (|^) prefers a negative a, 
which implies that one or more gaugino mass parameters could be negative (see Eq. @). We 
stress though that we do get some solutions with positive values of a as well. 

Now assuming, as an illustrative example, that A233 and A233 are the only dominant trilinear 
couplings, the factor K turns out to be X ~ Nc'm'^/Sir'^m-^. Taking the squarks and gauginos to 
be approximately at the 300 GeV scale, one obtains a ~ 2. 10-^ GeV-\ and ~ 1. lO'^ GeV. 
The minimum and maximum values of the input parameters that pass this test turn out to be 

A^33 ~ A^33 : [-1.3 X 10-3 , 1.3x10-3]; 

/^2,3 (GeV) : [-5.0 x 10-^ , 1.0 x 10-^]. (5) 



These values are consistent with the existing constraints on the above parameters [15, P^ . 

To conclude: If R-parity violation has to explain all neutrino data, it is essential to have 
both trilinear and bilinear ^ terms in addition to having a sterile neutrino in the model. Then 
it is possible to generate maximally mixed z/^ and I'r with their absolute masses in the eV range 
and mass-squared difference in the milli-eV^ range, as required by the data. The quoted ranges 
of the I? parameters within which we obtain solutions are based on certain simple-minded but 
plausible approximations made for the ease of presentation. The bottom line of our analysis is 
that we provide an affirmative answer to the question we have asked in the title. 
Note added : While we were finishing this note, we became aware of a preliminary SK solar 



neutrino analysis update [24| disfavouring a pure fg ^ I's solar neutrino oscillation at 95% CL. 
First, this is only a 2-a result which is not enough to exclude a model. Second, it has been claimed 
that a rate + spectrum combined analysis in a full four-flavour scenario exhibits an allowed zone 



[25|. So the fate of the sterile state may not be that dwindling, and we must wait till the SNO 



experiment tests this option. 

Both authors thank the CERN Theory Division, where this work has been done, for its warm 
hospitality. They also thank H. Dreiner and A. Raychaudhuri for illuminating comments on the 
manuscript. GB also thanks LPT, Orsay, for its hospitality, and S. Choubey and S. Goswami for 
clarifying remarks on neutrino data. 



'^In principle, the allowed spectrum should have been displayed as a five-dimensional plot where all the four 
input parameters are varying. Mostly for the purpose of a simplified presentation, merely to point out that there 
indeed exist solutions satisfying Eq. (^), we plotted the spectrum in two dimensions in Figs. 2a and 2b. In each 
plot, the 'other three' parameters are also varying. 
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Figure 1: The type lib four-neutrino mass pattern. Interchange of the members in a given pair or the 
relative location of the two pairs may be allowed. 
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Figure 2: Mass spectrum as a function of the (a) trilinear and (b) bilinear R-parity-violating parameters. 



